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La presente invention a pour objet une raethode pour deformer 
graduelleraent des representations ou realisations, generes par simulation 
sequentielle, d'un naodele stochastique non necessairement gaussien d'un 
milieu heterogene, en se basant sur un algorithme de deformation graduelle de 
5 modeles stochastiques gaussiens. 

La methode selon 1' invention trouve des applications dans la 
modelisation de zones souterraines ou il s'agit de generer des representations 
montrant comment est distribuee une certaine grandeur physique dans une 
zone du sous-sol (la permeabilite z par exemple), qui soient compatibles au 
"fo nueux avec des donnees observees ou mesurees : donnees geologiques, 
enregistrements sismiques, mesures obtenues dans des puits notanmient 
mesures des variations au cours du temps de la pression et du debit de fluides 
issus d' un gisement etc. 

Etat de la technique 

Par la demande de brevet FR 98/09.018, on connait une methode pour 
deformer graduellement un modele stochastique (de type gaussien ou 
apparente) d'un milieu heterogene tel qu'une zone souterraine, contraint par 
une ensemble de parametres relatifs a la strucmre du milieu. Cette methode 
comporte le tirage d'un nombre p (p=2 par exemple) de realisations (ou 
representations) independantes du modele ou d'une partie au moins du modele 
choisi du milieu, parmi Tensemble de toutes les realisations possibles et une 
ou plusieurs etapes iteratives de deformation graduelle du modele en 
effectuant une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes entre elles puis des realisations composites 
successivement obtenues tfventuellement avec de nouveaux tirages, etc.. les 
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coefficients de cette combinaison etant tels que la somme de leurs canes est 
egale a 1 . 

Les modeles gaussiens ou apparentes convienBent pour modeliser des 
champs de grandeurs continues et, de ce fait, sont mal adaptes a la 

modelis al ioa de z ones tr a versees p a r d e s re s eaux d e fractures qu de s s ystemes 

de chenaux, par exemple. 

On connait par ailleurs les algorithmes dits de simulation sequentielle, 
les plus largeraent utilises parmi les algorithmes de simulation geostatistique. 
Bien qu'ils conviexment particulierement pour la simulation des modeles 
10 gaussiens, il n'impliquent pas en principe une limitation a ce type de modele. 

On forme une representation geostatistique d'une zone souterraine par 
exemple en la subdivisant par un reseau a N mailles et en determinant un 
vecteur aleatoire a N dimensions Z = (Zi, Z2,-... Zn) correspondant au mieux a 
des mesures ou observations obtenues sur la zone. Conmoie le montre par 

15 exemple Johnson, M.E. ; in « Multivariate Statistical Simulation » ; Wiley & 
Sons, New York, 1987, cette approche reduit le probleme de la creation d'un 
vecteur a N dimensions a une serie de N problemes mono-dimensionnels. Un 
tel vecteur aleatoire n'est ni necessairement multi-gaussien ni stationnaire. La 
simulation sequentielle de Z implique d'abord la definition d'un ordre selon 

20 lequel les N elements (Zi, Zo Zn) du vecteur Z sont generes Ttm apres 

r autre. Sans rien perdre de generalite, on suppose que les N elements de Z 
sont generes en sequence depuis Zi jusqu'aZN. Pour tirer une valeur de chaque 
element Zi, (i = 1, N), il faut realiser les operations suivantes : 

a) construire la distribution de Zi conditionnee par (Zi, Z3.. ZuO 



F,(z,) = /'(Z, <^//Z,,Z, Z,_,) ;et 

b) tirer une valeur de Zi de la distribution Fc(Zi). 

Dans la pratique de la geostatistique, la simulation sequentielle est 
utilisee frequemment pour engendrer des vecteurs multi-gaussiens et des 
5 vecteurs d'indicateur non-gaussiens. La tache principale de la simulation 
sequentielle est la determination des distributions conditionnelles Fc(Zi) (i = 1, 
N). Des algorithmes et logiciels pour es timer ces distributions sont decrits 
par exemple dans : 

- Deutsch, C.V. et al, « uSLLb (Ueostatistical ilottware Library) and User's 
10 Guide » ; Oxford University Press, New York, Oxford 1992. 

En ce qui conceme le tirage des valeurs a partir de la distribution Fc(Zi), 
il existe aussi un vaste ensemble d'algorithmes connus. 



On considere la methode de distribution inverse par laquelle on obtient 
une realisation de Zi: = F~\u:) , oil XX{ est tire d'une distribution uniforme 
15 entre 0 et 1. D s'ensuit qu'une realisation du vecteur Z correspond a une 
realisation du vecteur U dont les elements Ui, Uo, ....,Un, sont mutuellement 
independants et uniformement distribues entre 0 et 1 . 

On s'apergoit qu'une simulation sequentielle est une operation S qui 
transforme un vecteur uniforme U= (Ui, U2, ....,Un) en un vecteur structure Z 
20 = (Zu Z3,.... Zn) : 

Z = S(U) (1) 



Le probleme que represente la contrainte d'un vecteur Z a differents 
types de donnees peut etre resolu en contraignant les distributions 
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conditionnelles Fc(Zi) (i = 1, N) et/ou le vecteur imiforme U= (Ui, Un, 
Un). 

Des travaux recents sur ralgorithme sequentiel ont ete concentres sur 
r amelioration de restimation des distributions conditionnelles Fc(zi) par 
mtegrauon de donnees geologiques et de donnees sismiques. 



On pent citer par example le papier de Zhu, H. et al ; « Formatting and 
Integrating Soft Data : Stochastic Imaging via the Markov-Bayes Algori thm » 
in Soares, A., Ed. Geostatistics Troia 92, vol.1 : Kluwer Acad. Publ., 
Dordrecht, The Netherlands, p. 1-12, 1993. 

10 Cependant, cette approche ne pent etre etendue a T integration de 

dormees non lineaires comme les pressions issues des tests de puits et les 
historiques de production, sauf a imposer une severe linearisation. En outre, 
conoume toute combinaison de vecteurs uniformes U ne donne pas un vecteur 
uniforme, la methode de deformation graduelle d'un modele stochastique 

15 developpee dans le brevet deja cite, ne peut etre directement appliquee dans le 
cadre de la technique sequentielle rappelee ci-dessus. 

Aussi, la methode selon T invention permet de rendre compatible les 
deux approches i.e. d'etendre le formalisme developpe dans le brevet 
precedemment cite, a la deformation graduelle de realisations, generees par 
20 simulation sequentielle, d'un modele stochastique non necessairement 
gaussien. 

Definition de la methode 



La methode permet de deformer graduellement une representation ou 
realisation, generee par simulation sequentielle. d'un modele stochastique non 



necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afm de la contraindre a un ensemble de donnees 
collectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa sti^cture. 

— ^ Elle est caractcrioGc en cc qu'cllc Lumpuiie T application d'un 

algoritiime de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) ayant un nombre N de variables mutuellement independantes qui 
est relie a un vecteur uniforme (U) a N variables uniformes mutuellement 
independantes par la fonction de repartition gaussienne (G), de fagon a definir 

up e chamc de realisations u(L) du vecLeiir (U). et 1' utilisation de ces 

realisations u(t) pour generer des realisations z(t) de cette grandeur physique 
que Ton cale par rapport aux donnees (non lineaires). 

Suivant un premier mode de mise en oeuvre, on definit la chaine de 
realisations u(t) du vecteur uniforme (U) a partir d'une combinaison lineaire 
J 5 de realisations du vecteur gaussien (Y) dont les coefficients de combinaison 
sont tels que la somme de leurs carres est egale a I'unite. 

Suivant un autre mode de mise en oeuvre. on realise une deformation 
graduelle d'un nombre n- de parties du modele representatif du milieu 
heterogdne tout en preservant la continuite entre ces n parties du modele, par 
20 une subdivision du vecteur uniforme (U) en un nombre n de sous-vecteurs 
mutuellement Lndependants. 

D'autres caracteristiques et avantages de la metiiode selon I'invention. 
apparaltront a la lecmre de la description ci-apres d'un exemple non limitatif 
de realisation, en se referant aux dessins annexes oix : 
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- la Fig.l montre la couche mediane d'une realisation d'un modele de facies 
servant de reference, generee par simulation sequentielle d' in die a trices ; 

- la Fig. 2 montre la variation en fonction du temps de la pression obtenue 
dans un test de puits du modele de reference ; 

5 - les Fig.3A a 3E montrent respectivement cinq realisations initiaJes de la 
couche mediane d'une zone reservoir, contraintes seulement par les facies le 
long du puit ; 

- les Fig. 4A a 4E montrent respectivement pour ces cinq realisations, les 
caui ' bes des pr e ssious du puits dans le modele de reference cuuipaiees a 

10 celles correspondant aux modeles initiaux ; 

- les Fig.5A a 5E montrent respectivement cinq realisations de la couche 
mediane du modele de facies conditionnees aux facies le long du puits et 
calees par rapport a la courbe de pression obtenue par test de puits ; 

~ les Fig. 6A a 6E montrent respectivement pour les cinq realisations, les 
15 courbes des pressions au puits dans le modele de reference comparees a 
celles correspondant aux modeles cales ; 

- les Fig.7A a 7E montrent respectivement comment les fonctions objectif 
correspondant respectivement a ces cinq exemples, varient avec le nombre 
d' iterations ; 
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- les Fig. 8A a 8E montrent les deformations graduelles operees par un 
changement de coefficient d' anisotropic, sur un modele a trois facies genere 
par simulation sequentielle d'indicatrices; et 




- les Fig. 9A a 9E montrent les deformations graduelles locales d'un modele a 
trois facies, generees par simulation sequentielle d'indicatrices. 

Description detaillee de la methode 

On considere une zone d' etude que Ton subdivise par une grille a N 
5 mailles et Ton cherche a construire des realisations ou representations d'un 
modele stochastique d'une certaine grandeur physique z representant par 
exemple la penneabilite des formations dans la zone. Le modele recherche doit 
s'ajuster a des doimees obtenues par des mesures ou observations en un certain 
nombre de points, et notamment a des donnees non lineaires. 

Le calage d'un modele stochastique a des donnees non lineaires pent 
etre formulee comme un probleme d' optimisation. On designe par f *'^= (fi°^' , 
p^obs ^ p^obs fp^^^)^ le vecteur des donnees non lineaires observees ou mesurees 

dans le milieu etudie (la zone reservoir), et par f = (fi , f2, fs, , fp), le 

vecteur correspondant des reponses du modele stochastique de la penneabilite 
Z = (Zi, Zo,..., Zn). Le probleme de contraindre le modele stochastique Z par 
les observations consiste k generer une realisation z de Z qui reduit a une 
valeur assez faible une fonction objectif que Ton deftnit comme la sonmie des 
erreurs quadratiques ponderees des reponses du modele par rapport aux 
observations ou mesures dans la zone-reservoir, soit : 

1=1 

ou CDi represente le poids attribue a la reponse fi. Les fonctions fi (i=L 2, 
p) ainsi que la fonction objectif O sont des fonctions du vecteur Z. On se 
trouve devant un probleme d' optimisation de dimension N. 





8 



Pour etendre le formalisme developpe dans le brevet deja cite, a la 
deformation graduelle de realisations generees par simulation sequentielle non 
necessairement gaussienne, on va partir d'un vecteur gaussien a N variables Yi 
avec i = 1, 2, N, mutuellement independantes, de moyenne nulla et de 
5 variance egale a Tunit^, et on va definir N variables uniformes mutuellement 
independantes Ui, Uo, U3,...Un par : 

Ui = G(Yi) V i = 1, 2, N 

ou G designe la fonction de repartition gaussienne standardisee. 

Ceci pose, on va appliquer ralgorithme de deformation graduelle 

10 developpe dans un cadre gaussien, au vecteur gaussien Y = (Yi, Y2, , Yn), 

pour construire une chame continue de realisations du vecteur uniforme U = 

(Ui, U2, Utj). Etant donne deux realisations independantes y.^^t yb de Y, on 

defmit la chame de realisations u(t) du vecteur U que Ton obtient selon la 
relation : 

15 w(0 = G(y^ cos/ + .v^ sinr) (2) 

Pour chaque t, u(t) est une realisation du vecteur U. Ensuite, par 

echantillonnage de la distribution conditionnelle F^{Z:) (i=l, 2, N) en 

utilisant les elements du vecteur u(t), on obtient un vecteur z(t) qui est. pour 
chaque t, une realisation du vecteur aleatoire Z. En consequence, on pent 
20 ajuster le parametre t comme dans le cas gaussien de fa^on a caler z(t) a des 
donnees non lineaires. Cette procedure est iteree jusqu'a atteindre un calage 
satisfaisant. 




Calage d'un modele de facies a des donnees de pression obtenues par 
des tests de puits 

Pour illustrer rapplication de la mdthode d' optimisation stochastique 
ci-dessus ddfinie, on va chercher a caler un modele stochastique de reservoir a 

^ des domiees de prcssions obtenues pai - des tests dt; puits. La Oonstiaiction du 

modele de reservoir est derivee d'une formation petroliere reelle qui comporte 
trois types de facies : deux facies reservoir de bonne qualite (facies 1 et 2) et 
un type de facies reservoir de tres mauvaise qualite (facies 3). La table 1 
definit les propridtes petrophysiques des trois facies : 











cD(%) 


c,(lQ-'har-') 


facies 1 


10. 


10. 


10. 


17 


2.1857 


facies 2 


1. 


1. 


1. 


14 


2.0003 


facids 3 


0.1 


0.1 


0.001 


9 


1.8148 
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Pour representer la distribution specifique de facies de la formation 
petrolifere, on genere d'abord one realisation binaire pour representer le facies 
3 et son complement. Ensuite, dans la partie complementaire du facies 3, on 
genere une autre realisation binaire independante de la premiere pour 
15 representer les facies 1 et 2. La formation est discretisee par un maillage 
regulier de 60x59x15 blocs de taille 15mxl5mxl.5m. Un modele de 
variogramme exponentiel est utilise pour estimer les distributions 
conditionnelles. La direction d'anisotropie principale est diagonale par rapport 
au maillage. Les portees du variogramme du facies 3 selon les u*ois directions 
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d'anisotropie sont respectivement de 300m, 80m et 3m. Les portees du 
variogramme des facies 1 et 2 selon les trois directions d'anisotropie sont 
respectivement de 150m, 40m et 1.5m. Les proportions des facies 1, 2, 3 sont 
respectivement de 6%, 16% et 78%. 

—5 On a realise un test de puits au moycn d' un simulateur do test de puits a 

differences finies, de la maniere decrite par : 

Blanc, G. et al : « Building Geostatistical Models Constrained by Dynamic 
Data - A Posteriori Constraints » in SPE 35478, Proc. NPF/SPE European 3D 
Reservoir Modelling Conference, Stavenger, Norway, 1996. 

10 On voit sur la Fig.l la couche mediane d'une realisation qui sert de 

modele de reference pour cette validation. La section du puits qui a ete 
perforee traverse horizontalement la couche mediane du modele de reservoir le 
long de I'axe x, Le diametre du puits est de 7.85cm, la capacite du puits est 
nuUe et les facteurs parietaux (skin) des facies 1, 2 et 3 sont respectivement 0, 

15 3 et 50. Le test de puits synthetique dure 240 jours avec un debit constant de 5 
m3/jour de fa^on a investiguer presque la totalite du champ petrolifere. La 
Fig.2 montre la variation de pression en fonction du temps. 

Le but etait de construire des realisations du modele de facies contraint 
par les facies rencontre le long du puits et par la courbe de pression obtenue 

20 lors du test de puits. La fonction objectif est deFmie comme la somme des 
differences quadratiques entre les reponses en pression du modele de reference 
et les reponses en pression de la realisation. Du fait que le comportement 
dynamique du modele de reservoir est principalement controle par le contraste 
entre les facies reservoir de bonne et de mauvaise qualites, la realisation 

25 binaire utilisee pour generer les facies 1 et 2 a ete fixee d'abord et seule la 
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realisation binaire utilisee pour generer le facies 3 a ete deformee pour le 
calage des donnees de pression. 

Les reponses en pression resultant des tests de puits pour les cinq 
realisations des Fig.3A k 3E, sont differentes de celle du modele de rdfdrence, 

commfi If, mnntrent les Fig. 1A h. 4 E. Partant lespecUvemem de ces b 

realisations independantes, en utilisant la medidde de calage iteratif ci-dessus, 
on obtient apr^s quelques iterations, cinq realisations calees (Fig. 5A a 5E) 
pour lesquelles les courbes de pression correspondantes sont tout a fait en 
accord avec celles du modele de reference, comme le montrent les Fig.6A a 
_LQ 6E. 
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Deformation gradueUe par rapport aux parametres structuraux 

Dans de nombreux cas, on ne dispose pas de donnees suffisantes pom- 
en deduire les parametres structuraux du modele stochastique : moyenne, 
variance, fonction de covariance, etc. Ces parametres strucuu-aux sont souvent 
15 donnes en termes d'intervalles ou de distributions a priori. Si leurs valeurs sont 
fausses, il est vain de rechercher une realisation calee sur des donnees non 
lineaires. Aussi, il est essentiel pour les applications de pouvoir proceder a une 
deformation gradueUe d'une realisation avec modification simultanee de 
nombres aleatoires et des paramed-es structuraux. L'algorithme de simulation 
sequentielle defmie par Tequation (1) foumit cette possibilite de changer 
simultanement ou separement I'operateur smicturel S et le vecteur uniforme U. 
Les Fig.8A-8E montrent les deformations graduelles que Ton obtient pour une 
realisation fixe du vecteur uniforme U quand on change le coefficient 
d' anisotropic. 
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Deformation graduelle locale ou regionalisee 

Quand les observations sont reparties sur des zones differentes d'une 
formation etudiee, un calage qui utiliserait mie deformation globale risquerait 
d'etre Lnefficace car 1' amelioration de I'accord que Ton obtiendrait sxir une 
zone pomrait Ic dctcriorcr sur une autre zone. D est done bieu pre fei able 
d'appliquer des deformations graduelles zone par zone. On considere une 

subdivision du vecteur U en un certain nombre n de sous-vecteurs U", 

mutuellement independants ce qui rend possible de realiser 
individuellement leur deformation graduelle. En appliquant separement a 
"TO ubacuxi des — suus-vecteuis U\ — U% — tf^; — ralyuiithme de — d^fon 



:uxi ues — suus-vecieuis — cn, — u% — en — l aiguatume ue — aeioimauou 

graduelle, on obtient une fonction de dimension n, du vecteur uniforme U : 



G(K;cosr,. + r;sin/J_ 



ou yj et Yl pour tout i =1, 2, n, sont des sous-vecteurs gaussiens 
independants. Pour un ensemble de realisations de Y^ et Kj, on resout un 
15 probleme d'optiniisation de n parametres ti, t^, tn pour obtenir une 
realisation qui ameliore, sinon maintient, le calage des donnees. Cette 
procedure peut etre iteree jusqu'a atteindre un calage satisfaisant. 

En realisant des deformations locales graduelles, on ameliore de fa^on 
significative la vitesse de calage dans tous les cas ou les mesures ou 
20 observations sont reparties dans differentes zones du milieu. 



• 
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On voit bien Teffet de cette deformation locale graduelle sur le modele 
a trois fades, des Fig. 9A a 9E on seule la partie inferieure gauche delimitee 
est affectee. 

La methode selon 1' invention pent etre facilement generalisee a la 
5 deformation graaueiie d une representaaon ou realisation d un modele 
stochastique quelconque, car generer une realisation d'un tel modele 
stochastique se ramene toujours a la generation des nombres uniformes. 
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REVENDICATIONS 



. 1) Methode pour deformer graduellemefll line lepieseuuQuu ou 

5 realisation, generee par simulation sequentielle, d'un module stochastique non 
necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
coUectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa structure, caracterisee en ce qu elle compone 1 apphcaduu J uu 

0 algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) h N variables mutuellement inddpendantes qui est relie a un 
vecteur uniforme U a N variables uniformes mutuellement independantes par 
la fonction de repartition gaussienne (G), de fa^on a construire une chaine de 
realisations u(t) du vecteur U, et I'utilisation de ces realisations u(t) pour 

15 aenerer des realisations z(t) de cette grandeur physique, que I'on cale par 
rapport aux donnees. 

0 2) Methode selon la revendication 1, caracterisee en ce que Ton definit 
une chaine de realisations u(t) du vecteur (U) a partir d'une combinaison 
lineaire de realisations du vecteur gaussien (Y) dont les coefficients de 
20 combinaison sont tels que la somme de leurs carres est egale a 1' unite. 

o 3) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle du modele representatif du milieu 
heterogene. simultanement par rapport aux parametres stmcmraux et aux 
nombres aleatoires. 




^ 4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle separee d'un nombre n de parties du 
modele representatif du milieu heterogene tout en preservant la continuite 
entre ces n parties du modele, par une subdivision du vecteur uniforme en n 
5 sous-vecteurs mutuellement independants. 




Fig.SA Fig.SB Fig.SC Fig.8D Fig.SE 




